Lattice study on kaon pion scattering length 
in the 7 = 3/2 channel 



Chuan Miao^, Xining Du^, Guangwei Meng^ and Chuan Liu' 

■ '^School of Physics 

■ Peking University 

O : Beijing, 100871, P. R. China 

(N 

^. 

Abstract 

00 . 

CN . Using the tadpole improved Wilson quark action on small, coarse and anisotropic 

lattices, Ktt scattering length in the 1 = 3/2 channel is calculated within quenched 
^ ' approximation. The results are extrapolated towards the chiral and physical kaon 

OO . mass region. Finite volume and finite lattice spacing errors are also analyzed and 

\ a result in the infinite volume and continuum limit is obtained. Our result is com- 

' pared with the results obtained using Roy equations, Chiral Perturbation Theory, 

■ dispersion relations and the experimental data. 
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1 Introduction 



Low-energy Kn scattering experiment is an important experimental method 
in the study of interactions among mesons [1,2,3]. It is also a good testing 
ground for our understanding of the low-energy structure of Quantum Chro- 
modynamics (QCD). Chiral Perturbation Theory [4,5], Roy equations [6], dis- 
persion relations [7,8] and other theoretical methods [9] have been used in the 
study of low-energy Ktt scattering. However, if the hadrons being scattered 
involve strange quarks, as in the case of Kn scattering, predictions within Chi- 
ral perturbation theory usually suffer from sizable corrections due to SU{3) 
flavor symmetry breaking, as compared with the case of pion-pion scatter- 
ing. Lattice QCD is a genuine non-perturbative method which can handle 
hadron-hadron scattering at low-energies. Therefore, a lattice QCD calcula- 
tion will offer an important and independent check on these results. In fact, 
pion-pion [10,11,12,13,14,15,16,17,18], pion-nucleon [11] and kaon-nucleon [19] 
scattering have been studied using lattice QCD. In this letter, we report our 
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quenched lattice results on the Kt: scattering length in the J = 3/2 channel. 
This quantity is of phenomenological importance for Kti scattering since the 
value of the scattering length in the 7 = 3/2 channel is correlated with that in 
the 1 — 1/2 channel and the value of the latter is crucial for the determination 
of the possible k resonance in the 1 — 1/2 channel [8]. 

The lattice gauge action that is used in our study is the tadpole improved 
clover action on anisotropic lattices [20,21]. Using this gauge action, glue- 
ball and hadron spectra have been studied within quenched approximation 
[20,21,22,23,24,25,26]. In our previous works, configurations generated from 
this improved action have also been utilized to calculate the tttt scattering 
lengths in the 1 = 2 channel [13] and the KN scattering length in the 1=1 
channel [19]. The fermion action used in this study is the tadpole improved 
clover Wilson action on anisotropic lattices [27,28,13]: 
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where various coefficients in the fermion matrix Ai are given by: 



r]i = u/{2us) , r]o = ^/2 , a = 1/{2k) - l/{lKmax) , 

Pt = C5w^(l + 0/(K) , Ps^csw/{2ut) . (2) 

Among the parameters which appear in the fermion matrix, csw = 1 is the 
tree-level clover coefficient and u is the so-called bare velocity of light, which 
has to be tuned non-perturbatively using the single pion dispersion relations 
[28]. In the fermion matrix (1), the bare quark mass dependence is singled 
out into the parameter a and the matrix A remains unchanged if the bare 
quark mass is varied. This shifted structure of the matrix A4 can be utilized 
to solve for quark propagators at various values of valance quark mass mo 
(or equivalently k.) at the cost of solving only the lightest valance quark mass 
value at K = «;max, using the so-called Multi-mass Minimal Residual {M^R for 
short) algorithm [29,30,31]. This is particularly advantageous in a quenched 
calculation since one needs the results at various quark mass values to perform 
the chiral extrapolation. 

The basic procedure for the calculation of Kn scattering length is similar to 
that adopted in the tttt scattering length calculation [13]. 
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2 Numerical calculation of the scattering length 



In order to calculate the elastic scattering lengths for hadron-hadron scatter- 
ing on the lattice, or the scattering phase shifts in general, one uses Liischer's 
formula which relates the exact energy level of two hadron states in a finite 
box to the elastic scattering phase shift in the continuum. In the case of Kt: 
scattering at zero relative three momentum, this formula amounts to a rela- 
tion between the exact energy E^^^ of the Kt: system with vanishing relative 
momentum in a finite box of size L with isospin /, and the corresponding 
scattering length Oq^ in the continuum. This formula reads [32] : 
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+ 0(L-«) , (3) 



where Ci = —2.837297, C2 = 6.375183 are numerical constants and iik-k ~ 
mKm-jr/ {rriK + ^tt) is the reduced mass of the Kir system. In this letter, the 
Kn scattering length Cq^^^^ in the isospin 1 — 3/2 channel will be calculated. 

To measure the hadron mass values m^r, rriK, rup and to extract the energy 
shift SE'^i;', one constructs the correlation functions from the corresponding 
one meson and two meson operators in the appropriate symmetry channel. 
For example, the Kn operator in the 1 — 3/2 channel is given by: 

o'KTit)-KitHit + i) , (4) 



where K^lt) and TTo{t) are the zero momentum kaon and pion operators 
respectively. The Ktt correlation function in the 1 = 3/2 channel then reads: 

c^Ktit) = {o'KTHt)o^Tm . (5) 



Numerically, it is more advantageous to construct the ratio of the correlation 
functions defined above: 

n'=\t) = it) /{CK{t)c„{t)) . (6) 



In this case, one uses the linear fitting function: 

^7=3/2(^)T>>t»l^_^^g/2)^ , (7) 



(3/2) 

to determine the energy shift SEj^'^ ' . Kir correlation function, or equivalently, 
the ratio TZ^^^^'^{t) is constructed from quark propagators, which are obtained 
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using the Multi-mass Minimal Residue algorithm with wall sources. Periodic 
boundary condition is applied to all three spatial directions while in the tem- 
poral direction, Dirichlet boundary condition is utilized. 

Configurations used in this calculation are generated using the pure gauge 
action for 6"^40, 8^40 and 10^50 lattices with the gauge coupling (3 = 1.9, 
2.2, 2.4, 2.6 and 3.0. The spatial lattice spacing is roughly between O.lfm 
and 0.4fm while the spatial physical size of the lattice ranges from 0.7fm to 
4.0fm. For each set of parameters, several hundred decorrelated gauge field 
configurations are used to measure the fermionic quantities. Statistical errors 
are all analyzed using the usual jack-knife method. Single pion, kaon and rho 
mass values are obtained from the plateau of their corresponding effective mass 
plots. The fitting interval is automatically chosen by minimal per degree of 
freedom. Due to the usage of finer lattice spacing in the temporal direction, 
good plateau behavior was observed in these effective mass plots. Therefore, 
contaminations from excited states should be negligible. These mass values 
will be utilized in the chiral extrapolation. 

After obtaining the energy shifts 6E^^^^\ these values are substituted into 
Liischer's formula to solve for the scattering length Oq^'^^ for all possible hop- 
ping parameter pairs: {k,u,k,s), which corresponds to the up (down) and the 
strange bare quark mass values, respectively. This is done for lattices of all 
sizes being simulated and for all values of (3. From these results, attempts 
are made to perform an extrapolation towards the chiral, infinite volume and 
continuum limit. 



3 Extrapolations of the scattering length 

The chiral extrapolations of physical quantities involving both the up (down) 
and the strange quarks consists of two steps. In the first step, the bare strange 
quark mass, or equivalently the corresponding hopping parameter Kg, is kept 
fixed while the hopping parameter of the up/down quark, Ku, is brought to 
their physical value k^'*^^. The precise value of be obtained by 

inspecting the chiral behavior of the pseudo-scalar (pion) and the vector me- 
son (rho) mass values. In the second step, one fixes the up/down quark hop- 
ping parameter at its physical value obtained in the first step and extrapo- 
late/interpolate in the strange quark mass. The physical strange quark hop- 
ping parameter k^^''^^ is determined by demanding the mass of the kaon being 
exactly at its physical value (the so-called K-input). 

We now come to the chiral extrapolation of the scattering length. In the chiral 
limit, the Ktt scattering length in the 1 — 3/2 channel is given by the current 
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algebra result: 

(3/2) _ .o^ 



where Oq^^^'* is the Ktt scattering length in the 7 = 3/2 channel and /jr ~ 
93MeV is the pion decay constant. To perform the chiral extrapolation of the 
scattering length, it is more convenient to use the quantity F — OQ^^^m'j,/ Hxn, 
which in the chiral limit reads: [13] 

(3/2) 2 -r 2 

F = ^ = - — -§ ^ -2.728 , 9 



where the final numerical value is obtained by substituting in the experimen- 
tal values for rup ~ 770MeV and ~ 93MeV. One-loop corrections reduce 
the central value of the scattering length by about 28%, with an estimated 
error of 40% [4,5]. The factor F can be calculated on the lattice with good 
precision without the lattice calculation of meson decay constants. Since we 
have calculated the factor F for several different values of valance quark mass, 
we could make a chiral extrapolation and extract the corresponding results in 
the chiral limit. 

In chiral extrapolations, wc have adopted a quadratic functional form and the 
fitting range of the extrapolation is self- adjusted by the program to yield a 
minimal psr degree of freedom. First, we extrapolate the factor F in the 
up/down quark mass values. Second, the extrapolated values for the factor F 
are then extrapolate/interpolate in the strange quark mass. In Fig. 1, we have 
shown the extrapolation in the strange quark mass for one of our simulation 
points. The resulting factor F after the up and down quark mass extrapolation 
are plotted in Fig. 1 versus the strange quark mass parameter 1/{2k,s). The 
line represents the quadratic extrapolation/interpolation and the final result 
is depicted as a solid square at the physical strange quark mass. The fitting 
quality Q for this fit is also displayed. 

After the chiral extrapolation, we now turn to study the finite volume effects of 
the simulation. According to formula (3), the quantity F obtained from finite 
lattices differs from its infinite volume value by corrections of the form 
In Fig. 2, we have shown the infinite volume extrapolation for our simulation 
points Bit P — 3.0, 2.6, 2.4, 2.2 and 1.9. The extrapolated results for the factor 
F are then used for the continuum limit extrapolation. 

Finally, we can make an extrapolation towards the continuum limit by elim- 
inating the finite lattice spacing errors. Since we have used the tadpole im- 
proved clover Wilson action, all physical quantities differ from their continuum 
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Fig. 1. Chiral extrapolation for the quantity F = OQ^mj,/ hkn for our simulation 
results at /? = 2.4 on 8^40 lattices. In this plot, quantity F, with the up and down 
quark mass already extrapolated to the chiral limit, is plotted (open squares) as a 
function of strange quark mass parameter 1/(2ks). The line represents the quadratic 
interpolation/extrapolation to the physical strange quark mass, where the result is 
depicted as a solid square. 
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Fig. 2. Infinite volume extrapolation for the quantity F 



,(3/2), 



rni/ hkit obtained 



from our simulation results at /? = 3.0, 2.6, 2.4, 2.2 and 1.9. The straight line 
represents the linear extrapolation in ag/vQ. The extrapolated results are shown 
with solid squares at = 0. 



counterparts by terms that are proportional to a^. The physical value of 
for each value of (3 can be found from Ref. [21,26]. 

The result of the continuum extrapolation is shown in Fig. 3 where the re- 
sults from the chiral and infinite volume extrapolation discussed above are 
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our simulation results at /3 = 3.0, 2.6, 2.4, 2.2 and 1.9. The straight line represents 
the linear extrapolation in ttg/rQ. The extrapolated result is also shown as solid 
square with the corresponding error. Also shown arc the result from current algebra 
and dispersion relations [7]. They coincide numerically and the value is show as an 
open circle. The chiral perturbation theory result is shown as a solid circle with 
error bar. 

indicated as data points in the plot for all 5 values of P that have been sim- 
ulated. The straight line shows the extrapolation towards the = limit 
and the extrapolated result is also shown as a solid square together with the 
chiral result from Ref.[4,5] which is shown as the open (tree-level, or current 
algebra result) and filled (one-loop) circles. Result from dispersion relations 
from Ref. [7] happens to coincide with the current algebra result (open cir- 
cle) numerically. It is seen that these results are compatible within error bars. 
To summarize, we obtain from the linear extrapolation the following result 
for the quantity F = OQ^'^^'m^p/ hk-k = (—1-87 ± 0.55). If we substitute in the 
physical values, we obtain the Kt: scattering length in the 7 = 3/2 channel: 
a^^^^rriT^ — (—0.048 ±0.014), which is to be compared with the one-loop chiral 
result of (—0.05 ±0.02) and the trec-lcvcl result (and also the early dispersion 
relation result) —0.07. Our lattice study also indicates a value consistent with 
the chiral result that is smaller in magnitude than the experimental result. 
The experimental results for the Kn scattering length tt'o^'^^mT, show sizable 
variation [1,2,3]. It lies between —0.13 and —0.05, significantly larger in mag- 
nitude compared with chiral and early dispersion relation results. In a recent 
analysis [8], a value of Og^^^-'m^ = — 0.129±0.006 was obtained using dispersion 
relations incorporated with chiral perturbation theory to fit the experimental 
data. 



4 Conclusions 

In this paper, we have calculated kaon-pion scattering lengths in isospin / = 
3/2 channel using quenched lattice QCD. It is shown that such a calculation is 
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feasible using relatively small, coarse and anisotropic lattices. The calculation 
is done using the tadpole improved clover Wilson action on anisotropic lat- 
tices. Simulations are performed on lattices with various sizes, ranging from 
0.7fm to about 4fm and with five different values of lattice spacing. Quark 
propagators are measured with different valence quark mass values. These en- 
able us to explore the finite volume and the finite lattice spacing errors in 
a systematic fashion. The infinite volume extrapolation is made. The lattice 
result for the scattering length is extrapolated towards the chiral and contin- 
uum limit where a result consistent with Chiral Perturbation Theory and early 
dispersion relation calculation is found. Our result for the scattering length is 
smaller in magnitude when compared with the experimental data results. 

Acknowledgments 

This work is supported by the National Natural Science Foundation (NFS) of 
China under grant No. 90103006, No. 10235040 and supported by the Trans- 
century fund from Chinese Ministry of Education. C. Liu would like to thank 
Prof. H. Q. Zheng and Prof. S. L. Zhu for helpful discussions. 

References 

[I] M.J. Matison et al. Phys. Rev. D, 9:1872, 1974. 

[2] N.O. Johannesson and J.L. Petersen. Nud. Phys. B, 68:397, 1973. 

[3] A. Karabouraris and G. Shaw. J. Phys. G, 6:583, 1980. 

[4] V. Bernard, N. Kaiser, and U.-G. Meissner. Nud. Phys. B, 357:129, 1991. 

[5] V. Bernard, N. Kaiser, and U.-G. Meissner. Phys. Rev. D, 43:2757, 1991. 

[6] P. Buettiker, S. Descotes-Genon, and B. Moussallam. hep-ph/ 03 10283, 2003. 

[7] C.B. Lang. Nuovo Cimento, 41A:73, 1977. 

[8] H.Q. Zheng, Z.Y. Zhou, G.Y. Qin, Z.G. Xiao, J.J. Wang, and N. Wu. Nud. 
Phys. A, 733:235, 2004. 

[9] B. Ananthanarayan, P. Buettiker, and B. Moussallam. Eur. Phys. J. C, C22:133, 
2001. 

[10] R. Gupta, A. Patel, and S. Sharpe. Phys. Rev. D, 48:388, 1993. 

[II] M. Fukugita, Y. Kuramashi, H. Mine, M. Okawa, and A. Ukawa. Phys. Rev. 
D, 52:3003, 1995. 



8 



[12] S. Aoki et al. Nucl. Phys. (Proc. Suppl.) B, 83:241, 2000. 

[13] C. Liu, J. Zhang, Y. Chen, and J.P. Ma. Nucl. Phys. B, 624:360, 2002. 

[14] JLQCD Collaboration. Phys. Rev. D, 66:077501, 2002. 

[15] K.J. Juge. hep-lat/0309075, 2003. 

[16] N. Ishizuka and T. Yamazaki. hep-lat/0309168, 2003. 

[17] CP-PACS Collaboration. Phys. Rev. D, 67:014502, 2003. 

[18] CP-PACS Collaboration, hep-lat/0309155, 2003. 

[19] G. Meng, C. Miao, X. Du, and C. Liu. hep-lat/0309048, 2003. 

[20] C. Morningstar and M. Peardon. Phys. Rev. D, 56:4043, 1997. 

[21] C. Morningstar and M. Peardon. Phys. Rev. D, 60:034509, 1999. 

[22] C. Liu. Chinese Physics Letter, 18:187, 2001. 

[23] C. Liu. Communications in Theoretical Physics, 35:288, 2001. 

[24] C. Liu. In Proceedings of International Workshop on Nonperturbative Methods 
and Lattice QCD, page 57. World Scientific, Singapore, 2001. 

[25] C. Liu and J. P. Ma. In Proceedings of International Workshop on 
Nonperturbative Methods and Lattice QCD, page 65. World Scientific, 
Singapore, 2001. 

[26] C. Liu. Nucl. Phys. (Proc. Suppl.) B, 94:255, 2001. 

[27] T. R. Klassen. Nucl. Phys. (Proc. Suppl.) B, 73:918, 1999. 

[28] Junhua Zhang and C. Liu. Mod. Phys. Lett. A, 16:1841, 2001. 

[29] A. Prommer, S. Giisken, T. Lippert, B. Nockel, and K. Schilling. Int. J. Mod. 
Phys. C, 6:627, 1995. 

[30] U. Glaessner, S. Guesken, T. Lippert, G. Ritzenhoefer, K. Schilling, and 
A. Prommer. hep-lat/9605008. 

[31] B. Jegerlehner. hep-lat/9612014. 

[32] M. Liischer. Commun. Math. Phys., 105:153, 1986. 



9 



